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Abstract 

We review the results obtained by three successful fluid dynamo ex- 
periments and discuss what has been learnt from them about the effect 
of turbulence on the dynamo threshold and saturation. We then discuss 
several questions that are still open and propose experiments that could 
be performed to answer some of them. 



1 Fluid dynamos 

It is now believed that magnetic fields of planets and stars are generated by the 
motion of electrically conducting fluids through the dynamo process. This has 
been first proposed by Larmor (Larmor 1919) for the magnetic field of the sun. 
Assuming the existence of an initial perturbation of magnetic field, he observed 
that "internal motion induces an electric field acting on the moving matter: and 
if any conducting path around the solar axis happens to be open, an electric 
current will flow round it, which may in turn increase the inducing magnetic 
field. In this way it is possible for the internal cyclic motion to act after the 
manner of the cycle of a self-exciting dynamo, and maintain a permanent mag- 
netic field from insignificant beginnings, at the expense of some of the energy of 
the internal circulation" (for reviews of the subject, see Moffatt 1978, Zeldovich, 
Ruzmaikin and Sokoloff 1983, Roberts 1994). 

The minimum set of parameters involved in a fiuid dynamo consists of the 
size of the flow domain, L, the typical fluid velocity, V, the density, p, the 
kinematic viscosity, v, the magnetic permeability of vacuum, /io, and the fluid 
electrical conductivity, cr. For most astrophysical objects, the global rotation 
rate, fi, also plays an important role. Three independent dimensionless param- 
eters thus govern the problem. We can choose the magnetic Reynolds number, 
Rm, the magnetic Prandtl number, Pm, and the Rossby number Ro 
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For planets and stars as well as for all liquid metals in the laboratory, the mag- 
netic Prandtl number is very small, Pm < 10~^. Magnetic field self-generation 
can be obtained only for large enough values of Rm for which Joule dissipation 
can bo overcome (for most known fluid dynamos, the dynamo threshold Rmc 
is roughly in the range 10 — 100). Therefore, the kinetic Reynolds number. 
Re = VL/v = Rm/Pm, is very large and the flow is strongly turbulent. This is 
of course the case of planets and stars which involve huge values of Re but is 
also true for dynamo experiments with liquid metals for which Re > 10^. Direct 
numerical simulations are only possible for values of Pm orders of magnitude 
larger that the realistic ones for the sun, the Earth or laboratory experiments. 
First because it is not possible to handle a too large difference between the time 
scale of diffusion of the magnetic field and the one of momentum; second, a 
small Pm dynamo occurs for large Re and requires the resolution of the small 
spatial scales generated by turbulence. Strongly developed turbulence has also 
some cost for the experimentalist. Indeed, the power needed to drive a turbulent 
flow scales like P oc pL^V^ and we have 



This formula has simple consequences: first, taking liquid sodium (the liquid 
metal with the highest electric conductivity), /iocr w 10m~^s, p w lO^kgm"^, 
and with a typical length scale L w Im, we get P ~ Rm', thus a mechanical 
power larger than 100 kW is needed to reach a dynamo threshold of the order 
of 50. Second, it appears unlikely to ever operate experimental dynamos at Rm 
large compared with Rmc- Indeed, it costs 8 times more power to reach 2Rmc 
than to reach the dynamo threshold. In conclusion, most experimental dynamos 
should have the following characteristics: 

• they bifurcate from a strongly turbulent flow regime, 

• they operate in the vicinity of their bifurcation threshold. 

Although the values of Rm and Pm that can be obtained in laboratory ex- 
periments using liquid sodium are not too far from the ones of the Earth core, it 
would be very difficult to perform experiments with large Rm at Ro significantly 
smaller than unity whereas we have Ro « 10~^ for the Earth core. The compar- 
ison is of course also difiicult in the case of the sun: although Ro is of order one 
for the solar convection zone, Rm is more than six orders of magnitude larger 
than in any laboratory experiment. As said above, the situation is worse when 
direct numerical simulations are considered. We thus cannot claim that cosmic 
magnetic fields can be reproduced at the laboratory scale except if we can show 
that the dynamics of the magnetic field weakly depends on some dimensionless 
parameters. 
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Relevant dimcnsionlcss parameters for some planets and the sun (from 
Zeldovich et al. 1983) and for laboratory fluid dynamos. {Rm has 
been evaluated on the full size even in the case of scale separation). 
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As already mentioned, laboratory dynamos operate in the vicinity of the 
instability threshold but at very high values of the Reynolds number. This give 
rise to a very interesting example of instability that differs in many respects 
from usual hydrodynamical instabilities. The dynamo bifurcation occurs from 
a base state which is fully turbulent. This may play a role on various aspects of 
the dynamo process and it raises several questions. 

• What is the effect of turbulent velocity fluctuations on the dynamo onset? 
Will they favor or inhibit dynamo action? Can they change the nature of 
the bifurcation? 

• Above onset, at which amplitude does the magnetic field saturate? Does 
it display an anomalous scaling with respect to the distance to threshold 

due to turbulent fluctuations? 

• What are the statistical properties of the fluctuations of the magnetic 
field? 

We will discuss this problem in connection with existing laboratory dynamo 
experiments. The first ones, performed in Karlsruhe and Riga, have been de- 
signed by taking into account the mean flow alone. Turbulent fluctuations have 
been inhibited as much as possible by a proper choice of boundary conditions. 
On the contrary, the VKS experiment has been first motivated by the study of 
the possible effects of turbulence on the dynamo instability. 

The paper is organised as follows: in section 2, we recall the results of 
Karlsruhe and Riga experiments. Section 3 concerns the effect of small turbulent 
fluctuations on the dynamo threshold and emphasises the effect of the Reynolds 
number on the scaling law for the saturated magnetic field above threshold. 
Section 4 considers some questions related to dynamos generated by strongly 
turbulent flows, i. e. with less geometrical constraints than the Karlsruhe and 
Riga experiments. The results of the VKS experiment with counter-rotating 
impellers are presented in section 5. A possible dynamo mechanism for the VKS 
experiment together with some additional comments are presented in section 
6. Finally, some open questions and other possible dynamo experiments are 
discussed in section 7. 
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2 The Karlsruhe and Riga experiments 



The first homogeneous fluid dynamos have been operated in liquid sodium in 
Karlsruhe (Stieglitz and Miiller 2001) using a flow in an array of pipes set-up 
in order to mimic a spatially periodic flow proposed by G. O. Roberts (1972), 
and in Riga (Gailitis et al. 2001) using a Ponomarenko-type flow (Ponomarenko 
1973). We first recall the flow geometries and briefly review the results obtained 
by both groups. 




(a) (b) 

Figure 1: Schematics for the experiments from Karlsruhe (a) and Riga (b) which 
show how helical flow is forced by guiding the sodium through steel channels 
(from Stieghtz and Miiller 2001 and Gailitis et al. 2001). 



2.1 The Karlsruhe experiment 

The experiment in Karlsruhe, Germany, was motivated by a kinematic dynamo 
model developed by G.O. Roberts (Roberts 1972) who showed that various 
periodic flows can generate a magnetic field at large scale compared to the flow 
spatial periodicity. One of the cellular flows he considered is a periodic array of 
vortices with the same helicity. Flows with such topology drive an a-effect that 
can lead to dynamo action. This mechanism is quite efficient at self-generation 
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(in the sense of generating a magnetic field at a low magnetic Reynolds number 
based on the wavelength of the flow). 

A dynamo based on this mechanism was constructed and run with success 
in Karlsruhe. A sketch of the experiment is shown in figure [1^. The fiow 
is located in a cylindrical vessel of width 1.85 m and height H = 0.7 m. It 
contains 52 elementary cells placed on a square lattice. Each cell is made of two 
coaxial pipes: an helical baffle drives the helical flow in the outer cylindrical 
shell whereas the flow in the inner shell is axial. In two neighbouring cells, the 
velocities are opposite such that the helicity has the same sign in all the cells. 
Although the volume is finite instead of the infinite extension assumed by G. O. 
Roberts, the dynamo capability of the flow is not strongly affected in the limit of 
scale separation, i.e., when the size L of the full volume is large compared to the 
wavelength I of the flow (Busse et al. 1996). In this limit, the relevant magnetic 
Reynolds number involves the geometrical mean of the two scales as a length 
scale, and the geometrical mean of axial and azimuthal velocities as a velocity 
scale. However, it can be shown using simple arguments that it is not efficient 
to increase too much the scale separation if one wants to minimise the power 
needed to reach the dynamo threshold (Fauve and Petrelis 2003). The flow is 
driven by three electromagnetic pumps and the axial and azimuthal velocities are 
independently controlled. The liquid sodium temperature is maintained fixed by 
three steam-evaporation heat exchangers. Measurements of the magnetic field 
were made both locally with Hall-probes and globally using wire coils. Pressure 
drops in the pipe and local velocity measurements were also performed. 

When the flow rates are large enough, a magnetic field is generated by dy- 
namo action. The bifurcation is stationary and the magnetic field displays 
fluctuations caused by the small scale turbulent velocity field (see figure [2]). 
This generation comes at a cost in the power necessary to drive the flow and 
the pressure drop increases. 

Due to the Earth's magnetic field, the bifurcation is imperfect but both 
branches of the bifurcation can be reached by applying an external magnetic 
field, as displayed in flgure [31 Among others, the experimentalists performed 
careful studies of the dependence of the dynamo threshold on the axial and 
helical flow rates and on the the electrical conductivity that can be varied by 
changing the temperature. They also considered the effect of flow modulation 
on the dynamo threshold and studied the amplitude and the geometry of the 
magnetic field in the supercritical regime. 

2.2 The Riga experiment 

The experiment carried out by Gailitis et al. (2001) has been motivated by one of 
the simplest examples of a homogeneous dynamo found by Ponomarenko (1973). 
A conducting cylinder of radius i?, embedded in an infinite static medium of 
the same conductivity with which it is in perfect electrical contact, is in solid 
body rotation at angular velocity fJ, and in translation along its axis at speed 
V . In an unbounded domain, this helical motion generates a travelling wave 
magnetic field. This Hopf bifurcation occurs for a minimum critical magnetic 
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Figure 2: Time recording of one component of the magnetic field in the Karl- 
sruhe experiment. The amplitude of the magnetic field increases after each 
increasing step of the flow rate of liquid sodium. Small fluctuations are visible 
once the magnetic field has saturated at a constant mean value. (Figure from 
Stieglitz and Miiller 2002). 

Reynolds number Rmc = fJ-oO'R^/ {R^)"^ + = 17.7 for an optimum Rossby 
number Ro = V/{Ril) — 1.3. We note that the maximum dynamo capability 
of the flow (Rmc minimum) is obtained when the azimuthal and axial velocities 
are of the same order of magnitude (Ro 1). This trend is often observed with 
more complex flows for which the maximum dynamo capability is obtained when 
the poloidal and toroidal flow components are comparable. 

The experiment set up by the Riga group is sketched in figure[TjD. Their flow 
is driven by a single propellor, generating helical flow down a central cylindrical 
cavity. The return flow is in an annulus surrounding this central flow. The 
geometry of the apparatus as well as mean flow velocity profiles have been 
optimized in order to decrease the dynamo threshold. In particular, it has 
been found that adding an outer cylindrical region with liquid sodium at rest 
significantly decreases Rmc- This can be understood if the axial mean flow as 
well as the rotation rate of the azimuthal mean flow are nearly constant except 
in boundary layers close to the inner cylinder. Then, the induction equation 
being invariant under transformation to a rotating reference frame and under 
Galilean transformations, the presence of some electrical conductor at rest is 
essential as it is in the case of the Ponomarenko dynamo. The three cylindrical 
chambers are separated by thin stainless steel walls, which were wetted to allow 
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Figure 3: The magnetic field amplitude increases above the critical flow rate 
in the Karlsruhe experiment. Another branch of self-generation can be reached 
only by imposing an initial field. This other branch is disconnected from the 
main branch. The imperfection of the bifurcation has been ascribed to the 
Earth's field. (Figure from Stieglitz and Miiller 2002). 



currents to flow through them. 

Figure m displays the growth and saturation of a time periodic magnetic field 
at high enough rotation rate. The nature of the bifurcation as well as dynamo 
growth rates have been found in good agreement with kinematic theory (Gailitis 
et al. 2002) that predicts a Hopf bifurcation of convective nature at onset. Note 
that this bifurcation can be affected by the ramping time scale of the propellor 
rotation rate (Knobloch 2007). In addition, the Riga group has made detailed 
observations of the magnetic field saturation value and the power dissipation 
needed to drive the fiow. These measurements give indications of the effect of 
Lorentz forces in the flow in order to reach the saturated state. It has been 
found that one effect of the Lorentz force is to drive the liquid sodium in the 
outer cylinder in global rotation, thus decreasing the effective azimuthal velocity 
of the inner flow and therefore its dynamo capability. Dynamo generation does 
also correspond to an increase in the required mechanical power. However, a 
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puzzling result is displayed in figure [S] the amplitude of the magnetic field for 
supercritical rotation rates does not seem to show the universal V — Rmc 
law. In addition, the form of the law seems to depend on the location of the 
measurement point. This is to some extent due to the absence of temperature 
control in the Riga experiment. Variations in temperature modify the fluid 
parameters (electrical conductivity, viscosity and density) and this should be 
taken into account by plotting the results in dimensionless form (Fauve and 
Lathrop 2005). 




100 200 300 400 500 

t[s] 

Figure 4: Time recording of the magnetic field from the Riga experiment. The 
dashed line gives the value of the rotation rate. The amplitude of roughly 
sinusoidal oscillations (not visible with the resolution of the picture) increases 
and then saturates when the rotation rate is increased above threshold (about 
1850 rpm). The amplitude of saturation increases when the rotation rate is 
increased further. The dynamo switches off when the rotation rate is decreased 
below threshold (Figure from Gailitis et al. 2001). 



2.3 Lessons from the Karlsruhe and Riga experiments 

Although there were no doubts about self-generation of magnetic fields by 
Roberts' or Ponomarenko-type laminar flows, these experiments have displayed 
several interesting features: 
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Figure 5: Magnetic field amplitude as the rotation rate is raised above the 
critical rotation rate. (Figure from Gailitis et al. 2003) 

• the observed thresholds are in rather good agreement with theoretical 
predictions (Busse et al. 1996, Radler et al. 1998, Gailitis et al. 2002) 

made by considering only the laminar mean flow and neglecting the small- 
scale turbulent fluctuations that are present in both experiments. 

• The nature of the dynamo bifurcation, stationary for the Karlsruhe exper- 
iment or oscillatory (Hopf) in the Riga experiment, is also in agreement 
with laminar models. 

• On the contrary, the saturation level of the magnetic field, due to the back 
reaction of the Lorentz force on the flow, cannot be predicted with a lami- 
nar flow model and different scaling laws exist in the supercritical dynamo 
regime depending on the magnitude of the Reynolds number (Petrelis and 
Fauve 2001). 

• Although secondary instabilities generating large scale dynamics of the 
magnetic field (such as field reversals for instance) have not been observed 
in the Karlsruhe and Riga experiments, small scale turbulent fluctuations 
of the magnetic field are well developed. 

These observations raise the following questions: 
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• What is the effect of turbulence, or of the magnitude of the Reynolds 
number, on the dynamo threshold Rmc'- Is it possible to observe how Rmc 
depends on Pm for a dynamo generated by a strongly turbulent flow (by 
changing Pm in experiments with a given flow at different temperatures 
for instance)? 

• What is the mechanism responsible for magnetic field fluctuations in the 
vicinity of the dynamo threshold: an on-off intcrmittcncy effect (Sweet et 
al. 2001) or chaotic advcction of the mean magnetic field by the turbulent 
flow? 

• What is the mechanism for field reversals? Is it possible to observe them 
in laboratory experiments? 

3 Effect of turbulence on the dynamo threshold 
and saturation 

3.1 Effect of small scale turbulent fluctuations on the dy- 
namo threshold 

As said above, dynamo experiments involve high Reynolds number flows and 
turbulent velocity fluctuations. Using the Reynolds decomposition, we write 

V(r,t)=V(r)+v(r,t), (3) 

where V(r) is the mean flow and v(r,t) are the turbulent fluctuations. The 
over-bar stands for a temporal average in experiments. The induction equation 
then becomes 

_ = V X (V X B) + V X (v X B) + j/n^V^B, (4) 

where = (/iocr) is the magnetic diffusivity. The Karlsruhe and Riga exper- 
iments have been designed by considering only the mean flow V(r) in order to 
predict the value of the dynamo threshold. However, we observe that turbulent 
fluctuations v(r, t) act as a random multiplicative forcing in the induction equa- 
tion (j4|). It is well known, both from simple theoretical models (Stratonovich 
1963, Graham and Schenzle 1985, Liicke and Schank 1985) and from experi- 
ments on different instability problems (Kabashima et al. 1979, Residori et al. 
2001, Berthet et al. 2003, Petrelis and Aumaitre 2003, Petrelis et al. 2005), 
that multiplicative noise generally shifts the bifurcation threshold. In addition, 
it sometimes modifles the dynamics of the unstable modes in the vicinity of 
threshold. In particular, multiplicative random forcing may generate intermit- 
tent bursting in the vicinity of instability onset (John et al. 1999, Berthet et 
al. 2003, Aumaitre and Petrelis 2003). This type of behaviour has been also 
observed with deterministic chaos instead of noise. It has been understood in 
the framework of blowout bifurcations in dynamical system theory (Piatt et al. 



10 



1993), and has been observed in a numerical simulation of the MHD equations 
without external noise (Sweet et al. 2001). Note however that in these simula- 
tions, Pjn is of order one, and the flow is chaotic at the dynamo threshold but 
not fully turbulent. We will use the same terminology "on-off intermittency" 
both for deterministic dynamical systems and systems with a noisy bifurcation 
parameter. 

The effect of velocity fluctuations has been considered a long time ago in 
order to explain reversals of the magnetic field of the Earth as a result of sta- 
tistical fluctuations of cyclonic convective cells (Parker 1969). Similar models 
using a noisy a-effect have been developed (Hoyng 1993). We will not study 
this type of problems here. 

As said above, no significant shift in threshold with respect to computations 
taking into account only the mean flow have been observed in the Karlsruhe 
and Riga experiments. Bursting phenomena in the vicinity of threshold have 
not been reported either despite the existence of velocity fluctuations. 

An explanation for the small shift in threshold has been given by Petrelis 
(2002) and Fauve and Petrelis (2003) as follows: in the limit of small fluctu- 
ations, one can calculate the threshold shift using a perturbation expansion. 
Let V be the average velocity field at onset, and B the neutral mode of the 
instability. Let V^^^ be the flow leading to the neutral mode B'^°) when there 
is no velocity fluctuations. Our aim is to find how the dynamo threshold of the 
velocity field V*^°) is modified in the presence of small turbulent fluctuations. 
We write v = i5v where i5 is a small parameter that measures the intensity of the 
turbulent fluctuations so that the amplitude of v is of order one. The neutral 
mode is likely to be slightly modified by the fluctuations as well as the dynamo 
threshold. We expand B and V in powers of S 

B^') are the corrections at order i to the neutral mode due to the presence of 
the turbulent fluctuations. Ci are constants that express the shift in the dynamo 
threshold caused by turbulence. We emphasise that we study the modification 
of the dynamo threshold of a mean flow with prescribed geometry due to the 
presence of fluctuations. When one inputs these expressions in equation the 
zeroth order part can be written 

L B W = - V X (vW X B W) - V^Bto) ^ ^ (5) 

L being the linear operator in the l.h.s. member. This is the laminar dynamo 
problem. By hypothesis, the instability onset is the one without turbulent 
perturbation. At next order in S we get 

L B(i) = ci V X (vC) X B^)) + V X (v X B^")) . (6) 

We now introduce a scalar product {f\g) and calculate the adjoint of L. As 
^^(0) — ^+ g^igQ nonempty kernel. Let C be in this kernel. Then 

(C|LB(i)) = (i+C|B(i)) = (7) 
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and this solvability condition gives the first order correction in the threshold 



(CIVx(vxBW)) 
' (C|Vx (V(o) xB(o))) ' ^ ' 

We use a scalar product in which the average over the realisations of the per- 
turbation is made. In that case, the average over the realisations of (C|V x (v x B*^" 
is proportional to the average of v, the value of which is zero by hypothesis. 
Thus, the dynamo threshold is unchanged up to first order in S, ci = 0. This 
result is obvious in many simple cases. For instance if v is sinusoidal in time, 
the threshold shift cannot depend on the phase which implies that it is invariant 
if V — » —V. This is also true if v is a random noise with equal probabilities for 
the realisations v and —v. Note however that simple symmetry arguments do 
not apply for asymmetric fluctuations about the origin although the threshold 
shift vanishes to leading order if the fluctuations have zero mean. 

To calculate the next order correction, we write equation ((4|) at order two in 
6 and get 

L B(2) = c2 V X (v(°) X BW) + V X (v X B(i)) . (9) 

We then get the second order correction 

^ (CIVx(vxBW)) 
' (C|Vx (V(o) xB(o))) ' ^ ^ 

where B^-*-^ is solution of 

LB(i) = V X (v X B(°)) . (11) 

Here, there is no simple reason for the correction to be zero. Its computation 
requires the resolution of equation (jlip . In some simple cases, an analytical ex- 
pression for C2 can be calculated and both signs can be found, thus showing that 
fluctuations can in general increase or decrease the dynamo threshold (Petrelis 
and Fauve 2006). 

We have thus obtained that when the amplitude d of turbulent fluctuations 
is small, the modification of the dynamo threshold is at least quadratic in S. 
Consequently, we can understand why the thresholds measured in the Karlsruhe 
and Riga experiments are very close to the predictions using the mean fiow V 
and thus ignoring turbulent fluctuations (the order of magnitude of the level of 
turbulent fluctuations related to the mean flow is certainly less than 10% in these 
experiments). The problem is more complex for experiments with unconstrained 
flows for which large scale fluctuations can be of the same order as the mean 
flow. We will consider it in the section about "turbulent dynamos" . 



3.2 Scaling laws for magnetic energy density in the vicin- 
ity of threshold 

In order to describe the saturation of the magnetic field above the dynamo 
threshold, we need to take into account its back reaction on the velocity field. 
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We thus have to solve the induction and Navier-Stokes equations that we restrict 
to incompressible flows (V • V = 0), 

— = V X (V X B) + V^'B, (12) 

at noa 

/d B"^ \ 1 

— + (V • V)V = -V - + — + lyV^Y + (B • V)B. (13) 

ot \p 2^0/ /io/O 

The flow is created, either by moving solid boundaries or by a body force added 
to the Navier-Stokes equation. We do not consider global rotation (fi = 0). 
We have to develop equations p2l [T3)l close to the dynamo threshold in order 
to derive an amplitude equation for the growing magnetic field. If the dynamo 
bifurcation is found supercritical, this allows the calculation of the saturated 
mean magnetic energy density {B'^)/2/j,q, where (•) stands for average in both 
space and time in this section. 

Thus, even in the simplest configuration, the problem involves three dimen- 
sionless parameters. One can choose, the magnetic Reynolds number, Rm, the 
magnetic Prandtl number, P™, and the ratio of the mean magnetic to kinetic 
energy density, leading to the following form of law 

ffi = p(y2)/(i?„,P,„). (14) 

Mo 

In general, the analytic determination of / using weakly nonlinear perturbation 
theory in the vicinity of the dynamo threshold is tractable only in the unrealistic 
case Pm ^ 1 such that the dynamo bifurcates from a laminar flow {Re ^ 1). 
For Pm ^ 1, a lot of hydrodynamic bifurcations occur first and the flow becomes 
turbulent before the dynamo threshold that occurs for Re ^ 1. We will briefly 
recall the scaling laws that are expected in both limits (for more details, see 
Petrehs and Fauve 2001). 



3.2.1 Saturation in the low Re limit 

At small Re, the Lorentz force should be balanced by the viscous force as it is 
larger than the inertial one in (fT3|) . When the magnetic field bifurcates above a 
critical velocity amplitude Vc, it thus generates a velocity modification of order 
SV given by 

SV_ _B2 

If this bifurcation is supercritical, we expect saturation for SV of the order of 
the distance to criticality V — Vc- In the bifurcation analysis, this balance results 
from the solvability condition. We thus obtain 

« _ (16) 

Analytic calculations in the small Re limit have been performed both for 
the Roberts flow (Gilbert and Sulem 1990, Busse and Tilgner 2001) and Pono- 
marenko type flows (Nunez et al. 2001). In the case of the Roberts flow for 
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which there are two spatial scales, the flow periodicity I and the size L of the 
full flow volume, the equation for (i?^) in the limit L ^ I takes the form 

(B^) CX ^ {Rra - Rmc), (17) 

with R,n = fio<7\/UVLl where U and V are respectively the axial and azimuthal 
typical velocities. The occurrence of I instead of L results from the balance 
between the Stokes and the Lorentz forces that both imply the small scale. The 
largest contribution to the current density is indeed the one associated to the 
small scale magnetic field in the limit of scale separation. 

More surprisingly, the basic property of this "laminar scaling", i.e., pro- 
portional to ly, or B^ / fiQpV^ oc 1/Re in dimensionless form, subsists in the case 
of some rotating flows (Childress and Soward 1972, Soward 1974). Equation 
(fT6|) is known as the "weak field scaling" of the Earth dynamo (Roberts 1988). 
The same property has been also found for a smooth helical flow in the limit 
Re > R„i > 1 (Bassom and Gilbert 1997). 



3.2.2 Saturation in the large Re limit 

At large Re, balancing the Lorentz force with the incrtial one gives 

VcSV ^ .^g. 
L pfioL ' 

Then, using SV oc V — Vc, we obtain the large Re scaling 

(^') « —7^1^ Rmc{Rm - Rmc). (19) 

The dimensionless group {B'^) p,Q{aL)'^ / p is known as the Lundquist number. 

The form of (|19p can be easily found by dimensional arguments. For Pm ^ 1, 
the Ohmic dissipative scale, 1^ = LRm^'^, is much larger than the Kolmogorov 
scale, Ik — LRe~^/^. Thus, the magnetic fleld grows at scales much larger 
than Ix and does not depend on kinematic viscosity. Then, we can discard the 
dependence of / on in (|14p . In the vicinity of dynamo threshold, V is not a 
free parameter but is such that ijlquVL « Rmc- We obtain 

{B^)'^^j^f{Rm). (20) 

The scalings for large and small Re differ by a factor Rmc/Pm ~ 10^ for 
experiments using liquid sodium (P^ « 10~^). It may be instructive to replace 
by the turbulent viscosity, i^t oc VL (respectively i^t VI in the case of scale 
separation) in the laminar scaling (|16p . Using V ~ Rmc/ PoO'L, we have 

(B^) CC ^ [Rm - Rmc) CX ,\ {Rm - Rmc) ■ (21) 

We thus recover the turbulent scaling. However, the above analysis does not 
require any assumption about the expression of the turbulent viscosity and is 
thus clearer. 
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3.2.3 Magnetic energy density of the Karlsruhe and Riga dynamos 

The Lundquist number, {B"^) 11^(0 rY / p where r is the radius of each cyhnder 
of the Karlsruhe experiment, is plotted as a function of in figure [S] for 
two different ratios of the axial to azimuthal flows. The large scale magnetic 
field being the dominant component and slowly varying in space, its local value 
thus gives a reasonable estimate of the mean energy density. The choice of the 
small scale r instead of the full size of the flow in the Lundquist number results 
from equation (jl7p . Rm has been defined using the geometrical mean of the 
axial and helical velocities as velocity scale but the height of the cylindrical 
volume H = 0.7 m has been taken as length scale. The later choice is somewhat 
arbitrary but allows a simple comparison with the other fluid dynamos. 

Figure [6] shows that the dynamo threshold depends on the axial to azimuthal 
flow ratio although its leading order effect has been taken into account in the 
definition of the velocity scale in i?„ . This may result from the expulsion of the 
transverse magnetic field by the azimuthal flow. Another possible explanation 
comes from the geometry of the helical flow. Indeed, part of the helical flow 
contributes to the axial flow. Taking this component into account in the defl- 
nition of Rm could decrease the difference between critical magnetic Reynolds 
numbers. On the contrary, the slope of the bifurcation curve does not seem to 
be strongly affected by the axial to azimuthal flow ratio. 
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Figure 6: Dimensionless magnetic energy density {B"^) fit^a^r^ / p as a function of 
the magnetic Reynolds number Rm for the Karlsruhe experiment. The different 
symbols are associated to different ratios of azimuthal to axial velocity. Rm is 
varied through change in the temperature leading to a variation of the electrical 
conductivity (data from Miiller et al. 2004). 



Figure [7| displays the dynamo bifurcation using the same parameters for the 



15 



Riga experiment. The Lundquist number as well as Rm are defined using the 
radius R = 0.125 m of the inner cylinder. The two curves are related to probes 
at different locations. We observe that the imperfection ascribed to the Earth's 
magnetic field in the Karlsruhe experiment does not exist in the Riga experi- 
ment. A constant external field is not resonant with the neutral mode that is 
a travelling wave and thus does not lead to any imperfection. Other measure- 
ments (not shown here) made at different spatial locations (see figure [5]) have 
displayed smaller values of the magnetic field that is concentrated in the vicinity 
of the shear for a Ponomarenko-type dynamo. Multiple point measurements and 
averaging thus should be made in order to have a better evaluation of the mean 
magnetic energy density in the flow volume. Fortunately, the difference between 
the predictions HH) and is so large that rough order of magnitude estimates 
of (B^) are enough. Taking into account the qualitative nature of our analysis 
in section 3.2, we conclude that the large Re scaling is in agreement with the 
experimental observations whereas the "laminar scaling" predicts a field that is 
orders of magnitude too small. We thus note that Karlsruhe and Riga exper- 
iments display an interesting feature: turbulent fiuctuations can be neglected 
when computing the dynamo threshold whereas the high value of Re has a very 
strong effect on the amplitude of the saturated magnetic energy density above 
the dynamo threshold. 



15 



A 

CM 

CQ 

V 



10 



34 



38 42 
Rm 



46 



Figure 7: Dimensionless magnetic energy density {B^)^oa^R^/p as a function 
of the magnetic Reynolds number Rm for the Riga experiment. The different 
values for same Rm are associated to different probe positions (data courtesy of 
Franck Stefani). 
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4 Laminar versus turbulent dynamos 



In a fully developed turbulent flow in a compact domain, i.e., without any mean 
flow prescribed by strong geometrical constraints, large scale; fluctiiaticms are 
generally of similar order of magnitude as the mean flow. Then, wo cannot a 
priori expect that the dynamo computed as if it were generated by V(r) alone 
can describe the observations as in the limit of small fluctuations considered in 
section 3.1. We thus propose the following definition: we call "laminar" a fluid 
dynamo which generates a magnetic field that displays the same characteristics 
as the one compiited as if the mean flow V(r) were acting alone. We mean by 
characteristics, the geometry of the mean field lines, the nature of the bifurcation 
(stationary, Hopf, etc) and the related dynamics of the large scale field, the 
approximate value of Rmc, etc. Otherwise, the dynamo is called turbulent, 
which means that turbulent velocity fluctuations v(r, t) qualitatively modify 
the nature of the bifurcation obtained with V(r) alone. This of course does not 
mean that V(r) has no effect on the dynamo process. This would remain true 
even if V(r) alone has no dynamo capability. 

4.1 Turbulent fluctuations and the shift of the dynamo 
threshold 

It has been believed that studying the dynamo generated by the mean flow can 
often provide a leading order description of a fluid dynamo even in the presence 
of large fluctuations. As just said, this is far from being obvious, but in that 
natural question is related to the effect of turbulent fluctuations on 
threshold as considered in section 3.1 in the limit of small fluctuations. There 
exists no general answer to this question. Even in the limit of scale separation, 
it is known that the role of turbulent fluctuations may be twofold: on one hand, 
they decrease the effective electrical conductivity and thus inhibit the dynamo 
action generated by V(r) by increasing Joule dissipation; on the other hand, 
they may generate a large scale magnetic fleld through the "a-effect" (Krause 
and Radler 1980, Moffatt 1978) or higher order similar effects even if V(r) = 0. 
However, most natural flows do not involve any clear-cut scale separation, and 
it has been recently shown in several numerical simulations that the "a-effect" 
may be very weak even for flow configurations involving large helicity where it 
could be expected to drive an efficient large scale magnetic field (Cattaneo and 
Hughes 2006, Hughes and Cattaneo 2007). 

The effect of turbulent fluctuations on the dynamo threshold has been also 
studied with direct numerical simulation of the Taylor-Green flow, but with 
different outcomes: the observed threshold has been found to increase in the 
presence of large scale fluctuations (Laval et al. 2006) whereas it has been 
observed to correspond to some threshold corresponding to the mean flow alone 
(Ponty et al. 2005, 2007). It should be emphasised that these direct simulations, 
as well as the previous ones about the "a-effect" , are not performed with small 
enough values of Pm so that one must be careful when using them to understand 
laboratory experiments. 
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In order to try to clarify this issue, it is instructive to consider the evolution 
equation of the magnetic energy 



where j is the current density {noj w V x B in the MHD approximation). At 
the dynamo threshold, the first term on the right hand side, the amplification 
term, should be equal to the second one that corresponds to Ohmic dissipation. 
If we assume that this is achieved in the absence of fluctuations, we observe that 
turbulent fluctuations can increase the dynamo threshold through two different 
types of mechanism: the most obvious one results from an increase of Ohmic 
losses due to the generation of magnetic fleld at small scales by advection of 
field lines by turbulent eddies. The second one is a loss of the efficiency of field 
amplification. The amplification term depends on the relative orientation of 
V, B and j, thus fluctuations in direction of these vectors can decrease it. We 
briefly discuss these two phenomena. 

The increase of Ohmic dissipation due to turbulence primarily depends on 
the value of the magnetic Reynolds number related to the corresponding ve- 
locity fluctuations, ^ia<^lv.rms where iv,,,,, if' the rms value of velocity fluctu- 
ations V and I is their integral scale. If nocrlvrms ~ 1 or smaller, the fluc- 
tuations of the magnetic field due to turbulence have a spectrum of the form 
cx {p,oaBo)'^k~'^\v\'^ oc {fj,oaBo)'^e'^^^k~^^^^ where Bq is the large scale mag- 
netic field and e = V^/L is the energy flux per unit mass (Golitsyn 1960, Moffatt 
1961). This spectrum being steeper than k~^, currents related to the large scale 
magnetic field are dominant and Ohmic dissipation is hardly increased by tur- 
bulence. The dynamo threshold is only slightly shifted as explained in section 3 
and observed in Karlsruhe and Riga experiments for which nocrlvrms is of order 
one or less. 

The situation strongly differs in flows without geometrical constraints for 
which Kolmogorov-type turbulent fluctuations exist up to the size L of the fluid 
container and have a rms value of the same order as the mean flow. The mag- 
netic Reynolds number of turbulent fluctuations is thus of the same order as 
the one related to the mean flow, and Rm > 10 at threshold. Then there 
exists an inertial magnetic range, = LRm^'^ < I < L, in which fol- 
lows a different scaling law. Unfortunately, no theoretical prediction for |_Bp is 
available in the vicinity of the dynamo threshold. Note that predictions for spec- 
tra of magnetohydrodynamic turbulence such as, |Bp cx (p^o)^^'*(e^o)"^^^fc~^^^ 
(Iroshnikov 1964, Kraichnan 1965) followed by many other models based on 
Alfven wave turbulence (for a review, see Verma 2004), require an applied field 
value Bq larger than the one achieved in the vicinity of the dynamo threshold. 
Similarly, the Kolmogorov-type spectrum |Bp cx p^Qe^/^k~^/^ , recently pro- 
posed for small P„i turbulent dynamos at large Rm (Fauve and Petrelis 2007), 
is not expected close; to threshold. 

Dimensional analysis lead to 




X. 



(22) 



\B\^ 0^ptlo{V^)k-^f{Rm,kla) 



(23) 
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where / is an arbitrary function. In the vicinity of the dynamo threshold it is 
expected to become 

OC pflo{V^)k-'{Rm - Rmc)g{kla) , (24) 

where g is another arbitrary function. Let us assume the following form of the 
spectrum of the magnetic fluctuations: 

= A(fc/fc^)" iffc/fc^<l 

= A(fc/A;<,)-"/3 iffc/fc^>l, (25) 

where a is an exponent that depends on the regime under consideration, k„ = 
and A is obtained by / \B\'^dk = {B"^). For a > —1, one gets A oc 
{B'^)LRm-^/'^ and for a<-l,A(x {B'^)LRm^"/'^. 
Then using 

-{f) = - I \j\'dk<x^ [ k'\B\'dk, (26) 
cr a J J 

we estimate Ohmic dissipation 



-(j^) o, ^^Rnv"/^ ifa>-l, 

(B^) Rm^^ ifa = -l 
HqCtL'^ Log(i?m) ' 

oc S^Rmi^+^-yA if _ 3 < c < -1, 



a \ ' if a < -3. (27) 



The next step is to relate {B"^) to (B). Following Moffatt (1961), one can 
assume that {B) imposes the value of the spectrum at kj, so that {B)^ « 
A{kLlkaTkL. Then (B^) oc (B)2Em(9+3")/4 if a > -1 and (B^) oc {Bf 
otherwise. One gets 



i(j2) oc if a > -3, 

oc if a < -3, (28) 

(B)"^ /{HQaL"^) being the dissipation in the laminar regime. We thus obtain 
that the effect of turbulent fluctuations is to increase Ohmic dissipation by 
a factor Rm"^ if a = 1/3 according to Moffatt (1961) or by a factor Rm^ if 
a = — 1 according to Ruzmaikin and Shukurov (1982). To sum up, any spectrum 
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less steep than in the interval [kL,ka-] leads to an enhancement of Ohmic 
dissipation. 

Fluctuations can also increase the dynamo threshold in the absence of any 

turbulent cascade. A simple example has been provided by Pctrclis and Fauve 
(2006) by considering phase fluctuations of the G. O. Roberts flow 

(V (cos {ky + (f>) — cos {kz + ip))\ 
U sin [kz + V') (29) 
U sin {ky + (t)) J 

where tjj and (p are two functions that depend on time only, or 

(V {cos {ky + (f)) — cos {kz + ^p)) \ 
U sin {kz + {V/k)d^^(t>cos {ky + (jj) (30) 
U sin {ky + </>) + {V/k)dxtp cos {kz + V")/ 

where (p and ijj depend on x only. Both flows involve phase fluctuations, i.e., 
motions of the eddies. In the first case, this amounts to switch the origin of the 
flow in time. In the second case, spatial fluctuations of the eddies along the x- 
axis are considered, the additional terms in the second and third components of 
the velocity field are just ensuring the incompressibility of the fiow. In the limit 
of scale separation, and assuming that the phase gradients in time (respectively 
ill enough, it has been shown that these large scale fliictuations 
always increase the dynamo threshold of the G. O. Roberts fiow. The above 
examples are simple enough to compute how fluctuations increase the threshold 
of the neutral mode generated by the mean flow alone. They also show that 
even if this threshold is increased, different modes that cannot be amplified by 
the mean flow, can be generated by the time dependent flow. It has been indeed 
shown that when ijj and cfi are sinusoidal functions of time in the first flow given 
above, fast dynamo modes are generated (Galloway and Proctor 1992). 

In conclusion, no general statement about the effect of fluctuations on the 
dynamo threshold can be made. When fluctuations occur at small scale and are 
of small amplitude such that their magnetic Reynolds number is small, their 
effect on the neutral mode generated by the mean flow alone is small and can 
be computed perturbatively. The shift in threshold occur at second (or higher) 
order. On the contrary, for Kolmogorov type turbulence with large scale fluctua- 
tions of the same order as the mean flow, it is expected that Ohmic dissipation is 
increased and the efficiency of the amplification mechanism of the neutral mode 
generated by the mean fiow alone can be decreased. However, other modes may 
be selectively amplified due to the presence of fluctuations. These observations 
question the validity of dynamo models based only on the mean flow, thus ne- 
glecting the effect of large scale turbulent fluctuations. Experimental studies of 
the transport of a localized magnetic field by a turbulent flow have also shed 
light on the effect of fluctuations (Volk et al. 2006b). In particular, they show 
a loss of the magnetic fleld orientation in the transport process that may lower 
the efliciency of the fleld amplification as stated before. 
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4.2 On-off intermittency 

In addition to a shift of threshold, the nature of the bifurcation can be modified 
in the presence of fluctuations. The simplest example is provided by the phe- 
nomenon of on-off intermittency. On-off intermittent behaviour is a common 
feature for bifurcating systems subjected to multiplicative noise. Consider a 
system close to an instability threshold. Multiplicative noise induces fluctua- 
tions of the instantaneous departure from onset. If the fluctuations are large 
compared to the mean departure from onset, the average over a time interval 
T of the instantaneous departure from onset can be negative even for long du- 
ration T. Then the amplitude of the unstable mode tends to zero (off-phase) 
before the average departure from onset turns back to a positive value and the 
unstable mode reaches amplitudes where nonlinearities saturate its growth (on- 
phase). This is a simple mechanism that leads to the random succession of 
phases where the unstable system is either close to its formerly stable solution 
or reaches values controlled by non-linearities. 

Up to now this behaviour has not been displayed by experimental dynamos. 
Two key elements have been identified that can limit the ability for a dynamo 
to exhibit an on-off intermittent magnetic field. 

The fluctuations that drive the off-phases are long time fluctuations. For 
simple dynamical systems, it has been shown that the on-off intermittent be- 
haviour is controlled by the zero frequency component of the spectrum of the 
multiplicative noise (Aumaitre et al. 2005, Aumaitre et al. 2006). Indeed, on- 
off intermittency takes place when the ratio between the departure from onset 
and the zero frequency component of the noise is small. It seems reasonable 
to assume that for dynamos, fluctuations of the velocity fleld at very small fre- 
quencies are required. It is probable that the currently working dynamos do 
not have enough low frequency fluctuations to generate an on-off intermittent 
magnetic field. 

Another effect that can play a role is the impcrfectness of the bifurcation. 
It is important for the mechanism of on-off intermittency that when the system 
evolves towards the off-phase, the effect of the noise term vanishes. It the bifur- 
cation is imperfect, the deterministic solution docs not tend to a zero magnetic 
field solution but to the imperfect branch. Then the effect of multiplicative 
noise does not vanish. This generates fluctuations that can destroy the inter- 
mittent regime (Pctrclis and Aumaitre 2006). Note that the effect is the same 
as for additive noise that also leads to a disappearance of on-off intermittency 
(Piatt et al. 1994). In the dynamo context, even though the fluctuations re- 
main multiplicative, the impcrfectness of the bifurcation can be the source of 
the fluctuations that prevent the observation of on-off intermittency. 

A source of imperfectness for experimental dynamos is the ambient magnetic 
field of the Earth and of the possibly magnetised parts of the experiment, as for 
instance the propellors in the VKS experiment. It might be helpful to screen 
these ambient fields in order to observe on-off intermittent dynamos. 

To sum up, in order to favour on-off intermittency above the dynamo onset, it 
is important to lower the source of imperfectness of the bifurcation, for instance 
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by screening the ambient fields and it is important to have a velocity field that 
displays large fluctuations at low frequencies. Then, slightly above the dynamo 
onset, on-off intermittency could be possibly observed. 

4.3 Sodium experiments with weakly constrained flows 

Early experiments have been performed by Lehnert (1957). A swirling flow has 
been generated by rotating a disk in a cylindrical vessel containing 58 1 of liq- 
uid sodium. Although a self-generated dynamo regime has not been reached, 
induction measurements have been performed with an applied axial field. The 
generation of a toroidal (w-effect) and poloidal field components have been ob- 
served. Another early sodium experiment has been motivated by planetary 
dynamos for which tidal forces have been considered as a possible source of 
power driving the dynamo (Malkus 1968). Precession of a rotating cylinder 
filled with liquid sodium has been studied by Gans (1970). The rotation rate 
was increased up to 3600 rpm for a precession rate of 50 rpm. An external mag- 
netic field amplification has been reported but no self-generation. 

Flows in spherical geometry are currently investigated in two dynamo exper- 
iments in the USA. The Wisconsin experiment (Forest et al. 2002, Normbcrg et 
al. 2006, Spence et al. 2006) is operational: the spherical shell is motionless and 
the fluid is driven by two propellors. This set-up is to be linked with numerical 
work by Dudley and James (1989). No dynamo action has been observed so far 
in the Wisconsin flow. Experiments in a similar geometry but at a smaller scale 
have been performed in Maryland (Peffley et al. 2000). 

Other sodium experiments, motivated by the context of the geodynamo, 
involve system-wide rotation. There is a widespread belief that global rotation 
may assist generation by lowering the critical Reynolds number. A spherical 
Couette flow has been studied in Grenoble. No self-generation has been observed 
but the effect of the Lorentz force due to an external magnetic field has been 
studied (Nataf et al. 2006). A similar instrument is being build in Maryland. 
A sphere of 3 meter in diameter, i.e., much larger than the currently running 
experiments, is being build with another concentric spherical shell inside. Both 
spheres can be put into rotation to investigate the effect of global rotation. 

5 The VKS experiment 
5.1 von Karman swirling flows 

The von Karman (VK) class of flows consists of flows in which the entrainment 
is performed by coaxial discs (Zandbcrgcn and Dijkstra 1987). The fluid is 
enclosed in a cylindrical shell. When the discs are operated in counter-rotation, 
these flows display various qualities of interest for a potential dynamo: a strong 
differential rotation and some helicity which arc key ingredients for a closed loop 
induction by u) and a effects. The flow lacks any planar symmetry which is also 
necessary for the a effect. As shown by measurements of pressure fluctuations, 
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large vorticity concentrations are produced (Fauve et al. 1993, Abry et al. 
1994) which may also act in favour of the amplification of the magnetic field 
if the classical analogy between vorticity and magnetic field production is to 
be believed. The choice of VK flows was thus motivated by the hope that the 
above features will make possible the generation a magnetic field by a strongly 
turbulent flow. 

The induction properties of these flows have been performed at moderate 
values of Rm up to 7 in Lyon to evidence the various induction steps that may 
lead to a dynamo (Odier et al. 1998, 2000, Volk et al. 2006a). With in mind the 
fact that, in the Riga and Karlsruhe dynamos, the observations fitted very well 
the prediction from kinematic dynamo computations, an optimisation procedure 
of the VK flow geometry has been implemented in CEA-Saclay in order to 
lower the threshold of the dynamo instability and to investigate a potential 
influence of the high turbulence level on the bifurcation (Marie et al. 2003, 
Ravelet et al. 2005). This procedure underlined in particular the importance of 
the electromagnetic boundary conditions, showing that electrically conducting 
propellors or possibly fluid motion behind them can increase Rmc by a factor 2 
(Stefani et al. 2006). 

5.2 The VKS2 experiment 
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Figure 8: Sketch of the VKS experimental set-up. The inner and outer cylinders 
are made of copper (in grey) and the other part in stainless steel. The dimensions 
are given in millimeter on the left half of the drawing and normalised by R on 
the right half. The 3D Hall probe is located either at point PI in the mid-plane 
or P2. In both cases, the probe is nearly flush with the inner shell. (Figure 
from Monchaux et al. 2007). 



The VKS acronym stands for "von Karman Sodium" . The VKS2 experiment 
is an evolution of a first design, VKSl (Bourgoin et al. 2002, Pctrelis et al. 2003) 
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which did not show any dynamo action. The changes compared to the VKSl 
design are the implementation of a coohng system and changes in boundary 
conditions motivated in part by the optimisation procedure discussed above. A 
sketch of the set-up is displayed in figure [51 The VK flow is generated in the 
inner cylinder of radius 206 mm and length 524 mm by two counter rotating 
discs of radius 154 mm and 371 mm apart. The additions to this base flow are a 
domain of sodium at rest surrounding the VK flow, an annulus in the mid-plane 
and pure iron discs. The sodium at rest has been shown to lower the threshold in 
the kinematic simulations using the mean flow. The annulus has been observed 
to stabilise to some extent the shear layer generated by the counter rotation of 
the discs. These changes in geometry were not sufficient to develop the dynamo 
instability and the last change was to use pure iron discs in order to modify the 
magnetic boundary conditions. In addition to the boundary conditions it also 
de-couples to some extent the domain in between the discs from the two domains 
behind the discs, whose flows may not be favourable at least in the kinematic 
simulations. This last configuration enabled the observation of a dynamo field 
as shown in figure [51 As the rotation rate of the discs is increased from 10 to 22 
Hz, one observes at the PI location (see fig. [H]) the growth of a magnetic field: 
the azimuthal component acquires a nonzero average value of order 40 Gauss 
with relatively strong fluctuations. The two other components display small 
average values but fluctuate with rms values of order 5 gauss. Even though the 
fluctuation level is much higher than in the Karlsruhe or Riga experiment, we 
call this dynamo stationary in the sense that it is not displaying any kind of 
time-periodicity or reversal. 

The magnetic Reynolds number in the VKS2 experiment is defined as Rm = 
Kfj^oaR^il where R is the radius of the cylinder and the rotation rate of the 
discs. K = 0.6 is a numerical coeflticient relating flR to the maximum velocity 
in the flow. With this definition the critical magnetic Reynolds number is close 
to 31 as can be seen from figure 1101 the circles correspond to the Lundquist 
number for VKS2 multiplied by a factor 25. 

5.3 Magnetic energy density of fluid dynamos 

The dimensionless magnetic energy of the three working dynamos is plotted in 
figure [TUI as a function of the magnetic Reynolds number. The striking feature 
is that all three dynamos bifurcate at critical Rm in between 30 and 35. Even 
though this may be pure coincidence due to a choice of definitions for it 
nevertheless means that critical values are really of the same order of magnitude. 

The saturation values of the three dynamos look very similar but again 
care has to be taken here as the Lundquist number is computed by using only 
a few localised measurements and not by integrating over the whole volume. 
This may cause the magnetic energy to be relatively badly estimated depending 
on the measurement point, even if one expects the scaling properties to be 
conserved. For instance, the ad-hoc factor C = 25 used for the VKS dynamo 
can be justified because measurements are performed at the boundary of the 
flow and thus underestimate the magnetic field intensity. 
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Figure 9: Time recording at position PI (fig. [5]) of the components of the 
magnetic field when the rotation frequency f2/27r is increased as displayed by 
the ramp below {R„i increases from 19 to 40). (Figure from Monchaux et al. 
2007). 



In any case the laminar scaling would predict magnetic energy 10^ too small 
such that for all three dynamos, the observations are compatible with the tur- 
bulent scaling of the saturation of the magnetic energy. 

For Rm > Rmc, the increase of the VKS2 Lundquist number seems of higher 
order than linear whereas the two other dynamos show an increase compatible 
with a linear trend. Whether this behaviour is reminiscent of some anomalous 
scaling is a difficult question to answer. The fact that the flow displays a much 
higher turbulence level than the two other dynamos may be a justification for 
anomalous scaling (see discussion below) but some care has to be taken in the 
interpretation of the data. In particular, the presence of the iron discs renders 
the bifurcation imperfect because of the remanence of the magnetisation in the 
discs (see discussion below). The apparent nonlinear scaling may be only due 
to the rather small range of available values for Rm — Rmc which may not allow 
to see the actual scaling. Also because of the imperfection, the apparent Rmc 
may appear lower than the actual one. 
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Figure 10: Lundquist number as a function of i?„j: Karlsruhe experiment for 
two ratios of the flow rates (■), Riga experiment (★) and VKS experiment for 
various working rotation rates and temperatures (•). C is a scaling factor used 
to display all curves on a single graph. C is one except for the VKS experiment 
for which C = 25. We do not expect to need this factor if measurements are 
performed in the bulk of the flow instead at the periphery as in Monchaux et 
al. 2007. 

6 Further comments and questions about the 
VKS experiment 

6.1 The effect of iron disks 

The addition of the iron discs enabled the growth of the dynamo field. Their ef- 
fect if first to modify the magnetic boundary conditions. For example, when the 
iron is still in the linear regime with permeability /io/i,- (i.e., the magnetisation 
is less that its saturated value), the ferromagnetic metal acts as a shield that 
prevents the field lines to go across the discs. Field lines are refracted when they 
penetrate the disc: the normal component of the magnetic field is continuous at 
the surface but the tangential one is increased by a factor /i^ in the discs. This 
may de-couple the main flow region in between the discs from the motion of the 
sodium behind the discs which has been shown by kinematic simulations based 
on the time averaged flow to increase the threshold of the kinematic dynamo. 

The magnetisation vector is also likely to be parallel to the plane of the discs 
in the same way as an elongated rod is more easily magnetised along its axis. 
One possibility that agrees with the axial symmetry of the experiment is that 
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the field lines of magnetisation would be loops centred on the cylinder axis. This 
would act in favour of magnetic field lines with the shape of loops parallel to 
the discs at least in their vicinity. Such loops with the same orientation as the 
magnetisation would be stabilised by the presence of the magnetised disc. Next 
to the disc, the flow is strongly outwards because of the centrifugal pumping. 
A material loop of small radius close to the disc would be stretched to a larger 
radius (see the sketch in figure[T2|). Thus a magnetic field loop of the same shape 
would be amplified by the stretching (Fauve and Petrelis 2003). The vicinity of 
the discs acts as an amplifier for the toroidal component of the magnetic field. 

As stated before, one consequence of the magnetisation of the discs is that 
the bifurcation is imperfect. The imperfection is of a different nature than the 
one in the Karlsruhe experiment. In the latter case, the imperfection was due 
to the Earth's magnetic field and it was possible to explore the second branch 
of solutions by imposing an additional external field as shown in figure [3l Here 
the situation is different: the magnetisation adds an additional contribution 
to the magnetic field but also the latter can modify the magnetisation. As 
iron has a low coercive force, the magnetisation can be flipped by a relatively 
small magnetic field. Let us illustrate this by a simple phenomenological model 
involving two coupled amplitude equations: 



The B variable is the analogue of the amplitude of the magnetic field and 
displays an instability for positive values of fi without the coupling to the second 
variable. M is analogous to the magnetisation. The solution A/ = is unstable 
even for i? = as the thermodynamically stable state is the saturated one and 
iron has a very narrow hysteretic cycle and low coercive force. Both variables 
are positively coupled as the presence of B induces some change of M due to the 
susceptibility and conversely any nonzero M causes a contribution to B. The 
nonzero stable solutions are displayed in figure[Tl] (for fj, > 1 other solutions exist 
but these are not relevant for the discussion here). One can see the imperfect 
bifurcation for B but also that both the positive and negative solutions display 
the same imperfection due to the coupling with the M variable. Indeed M and 
B always have the same sign because the low coercive force of iron causes the 
fiip of the magnetisation as soon as the magnetic field changes its sign. This 
very simple model illustrates the fact that contrary to the Karlsruhe dynamo, 
both branches of the imperfect bifurcation display the same feature. Applying 
an external field may force the dynamo to choose one specific branch but the 
shape would be the same in any case. 

Note that the B — 0, M = state is always linearly unstable. A more realis- 
tic description should incorporate space in order to describe magnetic domains 
and the metastability of the M = state. However, the aim of this model 
is just to illustrate that the sign of the imperfection depends on the one of 
the bifurcating magnetic field and thus is not prescribed by an external broken 
symmetry. 



dtB = nB-B^ + M, 
dtM = M- + B . 



(31) 
(32) 
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Figure 11: Stable solutions of the phenomenological model ([51]) . ([321) ■ On the 
left is the solution for B and on the right the one for M. Open and full symbols 
correspond to the same pair of solution {B, M). The dashed line correspond to 
the classical bifurcation without coupling B^ = /x. 



6.2 A possible dynamo mechanism for the VK flow 

Several predictions have been already presented for the VKS dynamo. They 
are all based on the mean flow alone; the neutral mode has been found dipolar 
with its axis in the mid-plane between the propellors, i.e. an equatorial dipole 
(Marie et al. 2003, Bourgoin et al. 2004, Ravelet et al. 2005, Stefani et al. 
2006). 

The geometry of the experimentally observed magnetic field differs from 
these predictions. We note that the mean flow being axisymmetric, it cannot 
generate an axisymmetric magnetic field according to Cowling's theorem. This 
constraint does not exist for the mean magnetic field generated by the full flow. 
We give below another possible mechanism that takes into account the helical 
structure of the radially expelled fluid within two neighbour blades which is 
averaged out when the mean flow is computed. 

The VK flows are very efhcient at converting the poloidal components into 
toroidal one by the u induction process: axial field lines are twisted by the 
strong differential rotation and this induces a toroidal component (Odier et al. 
1998). To get a loop back reaction from toroidal to poloidal, one can invoke 
an a effect localised close to the discs: the fluid ejected by the discs is strongly 
helical because of the shear induced by the differential rotation. This is the 
ingredient required for the a effect to occur. In this way, an axial field can be 
induced in the cylinder as drawn schematically in figure [T^ 

If one measures the magnetic field on the cylinder, the poloidal component 
is most likely weak as the currents that generate it are axisymmetric currents 
localised in the flow. On the contrary, the toroidal component should be stronger 
than the poloidal one and more likely to be even stronger when the measurement 
point is closer to the discs as observed in position P2. The observations are 
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Figure 12: A possible a — to mechanism for the VKS experiment using iron 
discs. Field loops are amplified by the outwards flow close to the discs. The 
fluid ejected by the discs is also helical due to the differential rotation and can 
convert the toroidal component Bg into a poloidal component Bp by a effect. 
The differential rotation converts also very efficiently the poloidal component 
into toroidal by ui effect. The relative signs of Bg and Bp depend on the sign 
of the helicity. 



compatible with this scenario but a deeper investigation of the dynamo field 
should be performed to validate it. 

In any case, an axisymmetric mean magnetic field cannot be generated by the 
VKS mean fiow because of Cowling's theorem. Thus, according to the definition 
given in section 4, the VKS experiment with counter-rotating disks generates a 
turbulent dynamo. 
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7 Open problems and other possible dynamo ex- 
periments 

7.1 Anomalous scaling laws for the magnetic energy den- 
sity? 

As said in paragraph 3.2, experimental dynamos operate with hquid metals for 
which Pm < 10~^ and neglecting this dimcnsionlcss number in the expression 
of the energy density leads to the large Re scaling (|20p . This seems reasonable 
because the magnetic field is dissipated below the Ohmic scale L/Rm^ which 
is much larger than the Kolmogorov scale L/Re^^'^ at small P„i- It is thus 
unlikely that the magnetic energy density depends on kinematic viscosity in 
this limit. However, one should keep in mind that there exist situations where 
a dimensionless parameter cannot be neglected even when it becomes much 
smaller (or larger) than the others. In these infrequent situations, the problem 
is said to be self-similar of the second kind if the dependence on this parameter 
is a power law (Barenblatt 1996). This assumption would give an additional 
dependence for the magnetic energy. 

We now want to discuss another "nonclassical" effect of the same class that 
may result from the presence of strong turbulent fluctuations of the velocity 
field. The magnetic field being forced by many different scales, one may expect 
a situation similar to critical phenomena for equilibrium phase transitions where 
large scale thermal fluctuations must be taken into account close to the critical 
point. This often leads to an "anomalous" behaviour of the order parameter 
in the vicinity of the transition, where it follows a power law with a critical 
exponent that differs from the mean field prediction. Although thermal fiuctua- 
tions may in principle also affect the amplitude of neutral modes in the vicinity 
of instability thresholds, it has been shown that this anomalous behaviour, if 
it exists, would be limited to a very small range and thus not detectable for 
most hydrodynamic instabilities (Hohenberg and Swift 1992). Various hydro- 
dynamical instabilities have been studied using liquid crystals (Rehbergei al. 
1991) or in the vicinity of the liquid-vapour critical point (Fauve et al. 1992, 
Oh and Ahlers 2003) in order to try to enhance the effect of thermal fluctua- 
tions. Although some effects have been observed, to the best of our knowledge, 
no anomalous dependence of the amplitude of unstable modes above threshold 
has been reported. Only mean fleld exponents, i.e., rational power laws given 
by simple symmetry arguments have been measured so far. For a supercritical 
instability, this assumption led to (jl9p . i.e., f{Rm, Pm) oc -y/ {Rm — Rmc), thus 
a mean field exponent 1/2. 

As said above, in phase transitions departure from mean-field theory oc- 
curs when the thermal fiuctuations cannot be neglected. This is expressed by 
the Ginzburg criterion that compares the amplitude of the order parameter 
predicted by mean field theory to the effect of thermal fluctuations (Ginzburg 
1960). Using a similar criterion, but taking into account the kinetic energy Ep 
of turbulent velocity fluctuations instead of /cT, would give an interval range 
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above the dynamo threshold where an anomalous scaling can be expected. It 
may be wrong to compare macroscopic turbulent fluctuations to temperature 
and estimating the number of relevant modes is difficult. However, if we assume 
that the mean field undergoes a pitchfork bifurcation and that the dependence 
of its amplitude in space involves a Laplacian to leading order, dimensional 
arguments lead to 



-°oso 



(33) 



where Bq is the pre-factor of the magnetic field above threshold in mean field 
theory and the one of its correlation length. Inserting the large Re scaling 
ifTO]) for Bq, and taking into account that fluctuations are of the same order of 
magnitude as the mean flow, we obtain 



Rm Rmc ^ 



(34) 



that can be of order one even if — L. We thus expect that a critical region in 
Rra with an anomalous scaling of the mean amplitude of the magnetic field can 
be easily observed. This can be understood since the energy of the fluctuations 
is signiflcant compared to the one needed to reverse the magnetic field slightly 
above the dynamo threshold. 

The above estimate of the range in — Kmc where an anomalous be- 
haviour can be expected depends of course on the nature of the bifurcation. 
For instance, if an equatorial mean field were generated, we expect the effect 
of fluctuations to be stronger than for an axisymmetric mean field because the 
orientation of the dipole can easily rotate under the infiuence of fiuctuations. 
This is reminiscent of the sensitivity of critical behaviour on the broken symme- 
try of the ordered phase: fiuctuations have a deeper effect when a continuous 
symmetry is broken. For instance the ordered phase is destroyed by the thermal 
fluctuations for a two dimensional XY spin model because the ordered phase 
breaks the axisymmetry of the problem. On the contrary, for a two dimensional 
Ising model, the ordered phase is axisymmetric and is not destroyed by the 
fluctuations (Goldenfeld 1992). 

We can argue that velocity fluctuations act in a multiplicative way whereas 
thermal fluctuations in phase transitions are additive. We note that in the 
presence of Earth's magnetic field, velocity fluctuations also involve additive 
forcing. In addition, even in the absence of any ambient field, the small scale 
fluctuations of the magnetic fleld are forced by an additive term resulting from 
the interaction between the small scale fluctuations of the velocity and the 
mean magnetic field generated by the dynamo. Differences from the mean field 
behaviour are expected as soon as the fluctuations of the magnetic field are of 
the same order as its mean value. Together with the hypothesis of self similarity 
of the second kind, the magnetic energy would then behave as 

( — ) = pV^P^,{Rm - R,ncf ■ (35) 

Mo 
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with /3 / 1/2. The results of the VKS experiment cannot presently be used to 
test this type of behaviour because of the imperfection due to the iron discs close 
to threshold. Thus, it is necessary to increase its maximum magnetic Reynolds 
number or to find a more efficient flow in order to obtain a turbulent dynamo 
without using iron discs. 

7.2 An optimised a-uj dynamo? 




Figure 13: An optimised a-io dynamo experiment. 

If one is to believe the induction scenario proposed above, then it may be 
possible to improve further the efficiency of the effects possibly involved in the 
VKS dynamo. We propose the set-up described in figure fT3l a VK flow is driven 
by two coaxial discs with a narrow gap. A very strong shear is thus created 
between these two impellers. The discs are fitted with blades resulting in a 
strong radial outward flow that expels the fluid. The discs centres are hollowed: 
the fluid expelled radially loops back into the gap by these openings. The radial 
velocity and the azimuthal shear result in a strongly helical flow. Therefore 
the flow displays two ingredients described above: a strong helicity and a strong 
differential rotation that drive respectively an a effect and an u) effect. With this 
modified set-up, one avoids the central region of the VKS experiment where the 
flow is radially inbound which tends to dampen the toroidal component of the 
flow. The expected unstable mode is qualitatively similar to the one described 
above for VKS, roughly axisymmetric and made of two parts: a poloidal field 
with same axis as the discs and an azimuthal field. The converging region is 
farther from the amplification region and is partially shielded if iron discs are 
used so that one expects this flow to be more efficient and to possibly lead to 
a lower critical magnetic Reynolds number. If the gain of efhciency due to the 
absence of radially inward flow in the mid-plane is large enough, one may also 
observe self-generation without using iron discs. 
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7.3 Non axisymmetric dynamos 

All currently working dynamos arc based on axisymmetric flows. The time- 
averaged velocity field of the Riga and VKS dynamos are invariant by rotation 
around the axis of their propellors. For the large scale magnetic field generated 
in the Karlsruhe experiment, the a-cfFcct created by the cellular flow is invariant 
by rotation around the axial direction of the flow. 

Symmetries usually do not favour dynamo action. There is no a-effect if the 
flow displays any planar symmetry. Anti-dynamo theorems are related to the 
symmetries of the flow or of the magnetic field. 

A non-axisymmetric flow can be generated by driving a fluid with a set of pro- 
pellors located at positions not aligned along their axis of rotation. A gallium- 
experiment based on such flow is currently studied at ENS-Paris (Berhanu, 
Mordant and Fauve 2007). It is driven by two pairs of propellors. Depending 
on their pitch and on the sign of their angular velocity, different flows that re- 
sult from interacting vortices generated by the propellors are observed. Such 
flows relieve the constraint of axisymmetry and may be efficient for induction 
or dynamo action. 

In addition, the large scale vortices interact, thus leading to fluctuations of 
the flow at low frequencies. This is probably important in order to generate 
dynamos that arc strongly affected by the flow fluctuations. Low frequency 
fluctuations are indeed required in order to have on-off intermittency in low 
dimensional dynamical systems (Aumaitre et al. 2005). Similarly, these fluc- 
tuations might favour the occurrence of anomalous exponents for the magnetic 
energy of turbulent dynamos. 

Finally we point out that a set of vortices generally leads to chaotic La- 
grangian trajectories. This is an important ingredient for fast dynamo action 
(Childress and Gilbert 1995). Although there is no hope to reach very large 
values of Rm in an experimental dynamo, studying the growth rate of the mag- 
netic field at moderate Rm in a flow displaying Lagrangian chaos at large scale 
can provide useful information. 

8 Concluding remarks 

The dynamo threshold and the geometry of the generated magnetic field in the 
Karlsruhe and Riga experiments were correctly predicted by taking into account 
the mean flow alone. This observation is supported by an argument showing 
that due to the low level of fluctuations, the shift in threshold has to be of 
order two at least in the amplitude of turbulent fluctuations. On the contrary, 
the geometry of the magnetic field generated in the VKS experiment cannot 
be generated by the mean flow alone. Large scale turbulent fluctuations are of 
the same order of magnitude as the mean flow because the flow is much less 
constrained. Although this needs to be confirmed with a detailed experimental 
study, we suggest a dynamo mechanism based on a a — w process that implies 
a non stationary part of the VK flow, suggesting that the observed dynamo is 
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turbulent (in the meaning used in this article). 

The three fluid dynamo experiments show that the magnetic energy density 
in the saturated regime above threshold is such that the Lundquist number is 
of the order of the distance to criticality. This results from the large value of 
the kinetic Reynolds number. 

The Karlsruhe and Riga experiments as well as the VKS experiment with 
propcUors in exact counter-rotation, have not displayed secondary bifurcations 
generating large scale dynamics of the magnetic field. However, it should be 
noted that for propellors driven at different rotation frequencies, a great variety 
of dynamical regimes (stationary, oscillatory, intermittent) as well as reversals 
of the magnetic field have been observed in the VKS experiment (Berhanu et 
al. 2007). It is likely that this results from the presence of competing instabil- 
ity modes that can be also observed in less turbulent dynamos by tuning two 
different control parameters. The effect of the fully turbulent character of the 
VKS2 experiment on these dynamical regimes deserves new experiments. 

We acknowledge Dr. F. Stefani and Dr. Robert Stieglitz for providing data 
and figures related to the Riga and Karlsruhe experiments (figures 1-7, 10). We 
thank all the participants of the VKS collaboration with whom experiments 
related to figures 8-10 have been performed. 
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